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Abstract. The two-component model of high temperature superconductors in its real 
space version has been solved using Bogoliubov-de Gennes equations. The disorder in 
q ' the electron and boson subsystem has been taken into account. It strongly modifies 

O . the superconducting properties and leads to local variations of the gap parameter and 

density of states. The assumption that the impurities mainly modify boson energies 
offers natural explanation of the puzzling positive correlation between the positions 
•^ of impurities and the values of the order parameter found in the scanning tunnelling 

iy~, ' microscopy experiments. 
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1. Introduction 

Recent years have witnessed the interesting series of experiments on high temperature 
superconductors (HTS). These include inter alia scanning tunnelling spectroscopy 
measurements which revealed atomic scale inhomogeneities on the surface of BiSrCaCuO 
PQ El El El E] and the quantum oscillations [SI E] which are signatures of the well 
developed Fermi surface. All these discoveries clearly indicate that after more than 
20 years of research, high temperature superconductors (HTS) still present a challenge 
and are the objects of intensive experimental [101 HH [12] and theoretical studies 

[HI El [IS1[I6]. 

The scanning tunnelling microscopy (STM) experiments are mainly performed 
on Bi family of HTS [17]. These materials cleave easily and it is straightforward to 
obtain atomically clean surfaces predominantly with BiO surface layer. This makes 
Bi superconductors also interesting materials to study by means of angle-resolved 
photoemission spectroscopy (ARPES) [9J. The strength of STM technique consists 
in its high energetic and spatial resolution. With STM one is able to study changes of 
local properties within the atomic distances. The poor screening, short superconducting 
coherence length of cuprate superconductors and low-dimensionality of the electronic 
structure contribute to large spatial fluctuations in normal and superconducting state 
properties. 

STM measures the quasiparticle excitations, existing at a given energy. From the 
differential conductance spectra one deduces e.g. local values of superconducting gap. 
The study of HTS reveals a number of unexpected features in the STM spectra. Among 
those "universally" observed there was nearly a homogeneous structure of conductance 
spectra at low energies despite strong variations at higher energies and in particular, 
inferred gap values [3], the asymmetry (i.e. G(— V) ^ G(V)) of the differential 
conductance [18] G(V) = dl/dV and characteristic dip- hump structure [19] . 

Already in the early days of HTS, the measured STM spectra have shown high 
degree of disorder [201 EI] ■ The subsequent more precise experiments allowed for the 
analysis of the results not only in real, but also in reciprocal space. It is important to note 
that with the Fourier transformed STM data, the possibility appeared to reconstruct 
the energy gap A(k), which turned out to be in agreement [22] with the ARPES 
measurements. The study of the spatial maps resulted in discovery that the values of 
the superconducting gaps are positively correlated with positions of oxygen dopants thus 
supporting the earlier claims that electronic inhomogeneities and atomic disorder are 
interrelated. Since then inhomogeneities are thus commonly attributed to the presence 
of disorder [23]. 

The development [21] of the so called "lattice-tracking spectroscopy" has enabled 
investigation of various correlations between the local system properties even at 
significantly different temperatures. This resulted in finding that not only the positions 
of dopant oxygen atoms correlate [51 [25] with local values of the superconducting 
gap, but that also the latter property is positively correlated with local density of 
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states (LDOS) measured in the normal state at a temperature well above that of 
the superconducting transition. This correlation indicates that normal state contains 
information about the interactions responsible for superconducting instability of the 
system [26]. 

The purpose of this paper is to study theoretically local properties of HTS relevant 
to the STM experiments. We describe the superconductor by the two-component boson- 
fermion (BF) model [27] [28]. For the homogeneous system the Hamiltonian can be 
written as 

k,cr <? 

+ 7^ E 9ta $f**frt + 4 t £ UA) ■ (1) 

k,q 

The cg(cl) operators refer to single fermions and annihilation (creation) operators 
of bosons are denoted by b^bV). The d-wave character of the order parameter in 
superconducting state is taken into account by assuming for two dimensional lattice 
g% = g[cos(k x a) — cos(k y a)]. 

Coexisting local pairs (bosons) and electrons are described by the Hamiltonian (TJQ) 
which has been proposed on phenomenological basis [27] but can also be derived from 
generalized periodic Anderson model [28J or as a low energy limit of cluster states [29]. 
Its relevance to HTS and other unconventional superconductors has been demonstrated 
in numerous works of various groups ^Dl SH S21 S3 E31 ESI ESI EZl SHI ESI SOI SH S2] . In 
a broader context this phenomenological two-component pairing scenario plays a role in 
a great variety of different physical situations, (i) It applies to resonant pairing in cold 
atomic fermion gases via a Feshbach resonance [13] , (ii) it captures the salient features of 
Anderson's resonating valence bond scenario [29] and (iii) it permits to draw conclusions 
about the pairing mechanism, if given by local dynamical lattice fluctuations, as was 
the original motivation for proposing this model [27] . 

As in this work we are concerned with the description of inhomogeneities and 
correlations observed experimentally, we generalize the Hamiltonian (pQ) to the real space 
and make the parameters random numbers. It means that disorder is responsible for 
various kinds of inhomogeneities. In the absence of ab-intio calculations [29J of how the 
boson and fermion parameters behave in an impure system [H] we study a few scenarios 
of their possible changes. 

It should be stressed that we are not averaging over disorder as experiments measure 
local values for a given (fixed) configuration. It turns out that with the BF model a 
number of experimental correlations found in STM can be understood with reasonable 
assumptions about the role of disorder. Here we shall be mainly interested in the maps 
of local density of states, the shapes of differential conductance curves as the functions 
of energy (applied bias) as measured at different sites and the values of energy gaps. 
Preliminary studies of this model have been presented recently [J5]. In literature there 
exist numerous attempts to describe theoretically various aspects of the STM spectra 
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in HTS cuprates including inter alia inhomogeneities [351 S3 SSI SHI E01 ED] , asymmetry 
of the differential conductance [52] and dip- hump structures [53J . 

The organization of the rest of the paper is as follows. In section 2 we introduce the 
two-component model, discuss the approach and compare the calculations for a small 
homogeneous system with bulk one. The results for the model with disorder in the 
electron and boson subsystems are presented in section 3 whereas main conclusions in 
section 4. 

2. The model and approach 

We start with the two-component model in a real space where fermions interact via 
short range forces with the hard-core charged and localized bosons. This boson-fermion 
model [281 EI] is represented by the Hamiltonian 

n BF = J2 %4A* + £ iyr - v) 4a« + £ W - 2 ^) 8 & 

i,j,(T icr i 



^ 2 



h \ (ciiCjt - CitCji) + bi {c]^ - ct+cjf)] , (2) 



<•./ 



where i and j denote the lattice sites, c ia (cj jCr ) stand for creation (annihilation) operator 
of fermion at the site i with spin a. b\ and bi are the creation and annihilation operators 
of hard-core bosons at the site i. ft stands for the chemical potential of the system and 
tij are hopping integrals, g^ is the electron-boson scattering (charge exchange). y^ mv is 
the local value of the fermionic energy level (scattering potential due to impurities) and 
Ef denotes the position of the bosonic level. We shall discuss the systems on a square 
lattice with hopping parameters between the nearest neighbors to be denoted t, which 
we use as the energy unit. 

2.1. Bogoliubov - de Gennes equations 

To account for the short coherence length and d-wave symmetry of the superconducting 
order parameter we assume that the boson-fermion coupling g^ takes on non-zero value 
for the nearest neighbor sites < i,j > only and equals +g if j — i ± x and -g if j — i ± y. 
Application of the standard Hartree-Fock-Bogoliubov decoupling of the type 
b\cilCjt ~ (bl)cilCjt + frl(cilA't) l ea ds to H BF = H B + H F , where H B is the single 
site bosonic Hamiltonian 

H B = £ (Ef ~ 2/i) b% + £ (blxi + kx*) (3) 

i i 

with the parameter Xi — /.(fti/2) < Ci^Cj^ — c^Cj^ > depending on the fermionic 

<j> 
degrees of freedom. H F is the Hamiltonian describing the disordered fermionic 
subsystem [54] . 

H F = £ tAc ia + J2 (Vr - A*) 4^ (4) 

i,j,(T icy 
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hi 

The standard statistical approach allows for an exact solution of the boson Hamiltonian 
H B . For a given realization of disorder one finds the following equations for the local 
number of bosons < b\bi > and the local boson average < b\ > 

<Pi> 4_S^ tanh (_|_), (5) 

<k> ^-gtanh(^). (6) 



We have denoted e, = \ J ( ' 2 M ) + |x«| 2 an d Ay = #y < b\ >/2. It is Ay which 

couples two subsystems as follows from equation (JHJ) below. 

The fermion part has the standard BCS-structure and we diagonalize it by the 
Bogoliubov-Valatin transformation, which introduces the new quasiparticle operators 
[55] %i, %^ and their adjoint counterparts jj^, 7/^ 

a a = Yl M^u + ^*7?,tl » ( 7 ) 



/ 



£ *,t = Yl H^'t - ^*7u 



This yields the following Bogoliubov-de Gennes equations pJ5] |5 




u\ 



E( A* ? ,-t,-(^-/i)^ II »'■ /"M „? )' (8) 



->y / \ ' ; / \ V; 



with Ay = Ay + Aji. For a given value of the chemical potential /x one solves the 
system of equations (jSJ) self-consistently together with (JBJ) for the energies E l and the 
functions w', v\. For a system on a square lattice ofnxm sites the size of the matrix is 
(2nm) 2 . If the total carrier concentration n = nf + 2n& is fixed, one finds jjl from 

n = TV E<4^} + 7I E < *& > ( Q ) 



and using equations flME]) and ffTT]) . 

The solutions enable calculation of all parameters of interest. In particular the local 
density of states (LDOS) N(E, i) at the site i as a function of energy E is given by 

N(E, i) = £ [\vi\ 2 S(E - E l ) + |^| 2 5(E + E 1 )] . (10) 

1 

For the purpose of numerical calculations we shall replace the Dirac 5(E) distribution 
by the Lorentzian ^ E ^ 2 , with r\ being a small smearing parameter. The local fermion 
number density n/^ = rif(Ri) = ^2 a (cf (T Ci (T ) ls found to be 

n f ,i = J2H 2 f(E l ) + H\ 2 (l-f(E 1 ))], (11) 
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where f{E l ) = l/[exp(E l /kBT) + l] is the Fermi-Dirac distribution function. Parameter 
Xi is given by 

* = EE^ HW - /(^)) +«K/(^)] ■ (12) 

i <j> 
The above equations are general in a sense that they allow for the study of both 
s-wave and d-wave superconductors depending on the choice of coupling gry. For the 
d-wave symmetry we define the local value of the order parameter at the site i as the 
staggered sum [50J over the neighboring bonds j, 



A; = (-A i>i+y + A iii+X - Ai^y + Ai^ x )/A. (13) 

The correlation function Cfh between various parameters /, h (in fact, also 
depending on the points of our lattice /, and hi) is defined as [49] 

/ / \Ji J average) \"j '''average ) 

Cf, h (R) = ■ I J - ■ ^, (14) 




where f average and h average are their values averaged over all sites and i? = |z — j\ is the 
fixed distance between the sites i and j in the above summations. 

In the next subsection we first recall the main features of the homogeneous boson 
- fermion model ([T]) and compare the results for a small cluster with those for the bulk 
homogeneous system. 

2.2. Clean system 

We consider here a real space version of the BF model and solve the corresponding 
Bogoliubov - de Gennes equations first for the homogeneous system of finite size with 
the periodic boundary conditions and compare the results with those obtained for the 
bulk system for the same set of parameters. This serves as a check of the quality of 
calculations for a relatively small size of the system. 

Figure Q] shows the dependence of the zero temperature order parameters (b), \x\ 
and the corresponding energy gap A (upper- left panel), the critical temperature T c 
(upper-right panel) and their ratio 2A/T C (bottom left panel) on the boson energy level 
E B measured with respect to the chemical potential. In the upper right panel the solid 
curve represents the data for the bulk homogeneous system, whereas the crosses show 
the results for the finite cluster. These data obtained for d-wave superconductor by 
solving the Bogoliubov - de Gennes equations for the very small system of size 17 x 21 
sites and the periodic boundary conditions agree very well with the results obtained 
from the mean-field study of the two- dimensional square lattice in the thermodynamic 
limit |H]. 

Small deviations are visible for E B > t, when real space calculations give slightly 
larger values. These are due to the finite size effect. Let us recall that the spectrum of 
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Figure 1. (Colour online) The local order parameters < bi >, \xi\ and (A^) at T=0 in 
the units t (upper left panel), superconducting transition temperature (T c /t) (upper 
right panel) and their ratio 2A/T C (lower left panel) are plotted as the functions of 
the boson energy {E B — 2/i)/t for a small, homogeneous system of the size 17 x 21 
with g=Q.375t, fi=0.0t, V lfnp =0.0t. The lower right panel shows the dependence of 
various (coupled) order parameters on temperature for E B =2/i. The superconducting 
transition temperature (crosses) obtained here is compared with that calculated for 
the corresponding bulk system (dotted curve). 



a small system is discrete and usually highly degenerate. The degree of degeneracy, 
however, is substantially reduced if one chooses rectangular shape of the system. 
This reduces number of symmetries in the Bogolubov-de Gennes equations and lowers 
degeneracy of energies. Best agreement with the bulk data is obtained for lengths n and 
m being different prime numbers. 

We remark that the quantities plotted in Figure Q] are symmetric with respect to 
E B — 2/i = 0. Physically this property comes from the hard-core nature of the local 
pairs (LP). Upon filling the bosonic level we obtain a perfect particle- hole symmetry 
between n B and 1 — n B occupancies. If /x = both these cases equally promote 
the superfluid/superconducting abilities of the coupled boson-fermion system (J2J). In 
particular, let us notice strong enhancement of the critical temperature near the half- 
filled boson level which is accompanied by the non-BCS relation 2A/T C (marked in lower 
left panel of Fig. [TJas LP) being a hallmark of the present scenario [27J. 

The model shows superfluid characteristics which are intermediate between those 
of a BCS regime and those of local bipolaronic pairs [30j ES] . For the homogeneous case, 
the effective interaction 7q between fermions depends on the boson-fermion coupling 
g, the position of the bosonic energies E B with respect to the chemical potential and 
hard-core boson concentration ns as Iq = (g 2 /(E B — 2/x))(l — 2ub)- This can be 
easily seen (391 HI] by rewriting equations (JSJ) -(jHJ) for the clean bulk system in the form 
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of effective BCS equation for the fermionic order parameter. The coupling g^ induces 
BCS like pairing in the fermionic subsystem and the Bose-Einstein condensation (BEC) 
in the bosonic subsystem. The inherent property of the model is that superconductivity 
vanishes at the same temperature T c in both subsystems [12] . Temperature dependencies 
of the bosonic | < b > |, fermionic \xi\ and Aj order parameters of the clean system are 
shown in the lower right panel of Figure [U 

For the vanishing boson-fermion coupling g, t j = model (T5]) describes the two 
non-interacting subsystems; dirty fermions and immobile hard-core bosons in random 
potential. None of the subsystems taken alone undergoes the transition. 

2.3. Remarks on the phase diagram of HTS 

The undoped parent compounds of cuprate HTSs are characterized by a non- 
superconducting ground state, which is antifferromagnetic insulator. With the 
increasing doping the system starts to be metallic and superconducting. The 
superconducting transition temperature T c increases with increasing doping, attains 
a maximal value for the optimal doping and decreases again in the overdoped regime. 
At the same time the normal state properties change dramatically. The pseudogap state 
[57] , which sets in at small doping disappears and the overdoped system might behave as 
normal Fermi liquid. The doping is on the one hand the source of charge carriers in the 
model, while on the other hand it introduces impurities. Judging from the appearance 
of superconductivity for finite impurity concentration the plausible assumption would 
be that impurities modify local values of the parameters in the Hamiltonian fl2]) in such 
a way as to promote superconductivity. 

We thus assume that oxygen impurities in HTS are the source of carriers (holes) 
and disorder in the system. The previous studies of the disordered boson-fermion model 
concentrated on the bulk properties and required the averaging over configurations. 
This was achieved by averaging the free energy [58] or using the coherent potential 
approximation to calculate the averaged Green functions [541 159] . The role of disorder 
is hard to overestimate in these materials. It seems that most of their unusual properties 
may simply be due to an unprecedented degree of disorder resulting from both potential 
scattering and random values of other parameters. The simultaneous presence of strong 
correlations in real materials makes the problem more complicated [60] even at the mean 
field level [61]. 

In the two-component model it is the boson-fermion scattering which induces 
superconducting transition. In this scenario the pseudogap phase is dominated by 
the incoherent pairs of bosonic character. To model the initial increase of the 
superconducting transition temperature with doping we assume that dopants modulate 
bosonic levels Ef changing them locally from the value well below or above the 
Fermi level to its vicinity which allows for efficient fermion-boson scattering and 
superconducting instability. The proper choice of Ef leads to positive correlations 
between the positions of impurities and the amplitude of the gap in the local density of 
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states. 

3. Results of calculations 

The description of the cuprate superconductors by the disordered boson-fermion model 
(J2]) allows for a great variety of scenarios. In the following we shall consider a number 
of them as our main goal here is to look at various possibilities offered by the model at 
hand. 

We numerically solve the Bogoliubov-de Gennes equations for small clusters with 
the size n x m sites. We use periodic boundary conditions, with which the solutions 
converge to the correct bulk values, even for a relatively small system. We determine the 
values of the order parameter Ay for the d-wave superconductivity at each bond i — j, 
the carrier concentration n, = n(f*j) at each site and the local density of states N(E, i). 
The results are presented in the form of maps showing these local parameters. For the d- 
wave order parameter we shall use the local "representation" of A^- as given by equation 
( TT3l) . In this work we limit the discussion to the analysis of the ground state (T = OK) 
properties. In the following we shall compare our results with STM data taken at low 
temperatures. Because the energy gap changes very slowly for temperatures Tj0.3T c , as 
it is visible in Fig. 1 the ground state data can safely be compared with the experiment 

As mentioned before, the model and approach allow us to study various kinds of 
inhomogeneities. They may result from the local changes of: 

- Ef , the position of the boson energy level. As we know from the study of the 
model for the clean system, the position of Ef with respect to the chemical potential 
has a profound effect on the appearance of superconductivity. 

yimp ^ yjYiicYi can be interpreted as local atomic levels. This term represents the 
scattering of fermions by point impurities. 

- tij, the fluctuating values of hopping parameters. The fluctuations of this type 
are not discussed in this paper. 

3.1. Random bosonic levels supporting or suppressing superconductivity 

It is the experimental fact that undoped copper oxides are not superconducting. 
As noted above, the existence of coupling between the subsystems g^ causes 
superconducting instability. It is also evident from Figure [1] that for superconductivity 
to appear the bosonic levels Ef have to be close enough to the Fermi level. It seems thus 
natural to expect that doping of parent compounds introduces carriers into the system 
and places the boson levels in the close vicinity to the Fermi energy. In accordance with 
such a picture we study here a few different (model) scenarios of doping. 

First, we consider random Ef centers, which we also call impurities with parameter 
values supporting superconductivity. It means they move locally the boson level towards 
the Fermi level (c.f. upper panel of Figured]). Second case is another extreme in which 
contrary to the above expectations impurities move the boson level out of its initial 
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position in the close vicinity of \x and push it up or down on the energy scale. This 
makes the boson-fermion scattering ineffective due to phase space restrictions. In this 
scenario Ef impurities suppress superconductivity. 

The above two scenarios neglect the effect of disorder on the fermion subsystem i.e. 
we put V* mp = 0. In real systems one expects that the centers with the boson levels 
closer to the Fermi energy, at the same time are the source of (presumably strong) 
potential scattering for fermions, which in our system is modelled by V imp . Thus as the 
third scenario we consider the system in which the doping will introduce randomness 
in both boson and electron on-site parameters. In order not to increase a number of 
parameters we assume [62] that in this case the relation Ef = 2Vi and V imp = Vi is 
valid. 

It has to be stressed that in principle there exists an additional degree of freedom 
which is related to the spatial extent (50] at which introduction of the impurity into the 
lattice will modify the parameters of the effective model as d2J. Accordingly we discuss 
short (local) and long range (extended) impurities. In the latter case the impurity at a 
given site will modify the boson energy levels in neighboring and more distant sites (see 
Figs. EH6]). If two extended impurities happen to influence the same site their effects 
simply add up. 

In Figure [2] we show a superconducting system with relatively low transition 
temperature and introduce local impurities which move the boson levels Ef towards 
the chemical potential, thus increasing T c . Next two figures (i.e. Figure [3] and H]) refer 
to the extended impurities. The maps of the local order parameter Aj of the d-wave 
superconductor are shown in the upper left panels. 

The upper right panels of Figures [3] and 0] show the distance R dependence of the 
correlation functions Cf t h, as defined in Eq. (114J1 . for the following parameters (/, h): 

(i) local values of the fermionic order parameter Aj with positions of the impurities; 
R im P, 

(ii) local values of bosonic order parameter (6j) with R %m v 

(iii) Aj with the maximal height of the local density of states N(i) denoted as 
LDOS peafc in the figures. Note, that at each lattice site the height of N(i,u) is the 
largest at slightly different energy oj. 

(iv) the correlation of A.; with its value at the distance i + R. The extent of this 
correlation can be viewed as a measure of the coherence length £. For the short range 
impurities it is of the order of one lattice spacing and increases with the range of effective 
interaction. 

The STM spectra shown as the functions of energy (bias) in Fig. [2] are taken 
along vertical line at X = 53 (left) and X = 67 (right), whereas those in Fig. [3] are 
for X = 17 and X = 43, respectively. The size of the system is 71 x 77 and there are 
16% of impurities (marked by circles) which change the boson levels Ef locally from the 
homogeneous value 0.58t to Ot (Fig. [2]). Other parameters of this d-wave superconductor 
are t\ = l.Ot, g = 0.5t, V t mp = 0, and the total number of carriers n = 1.30. 

It has to be noted that even without impurities the homogeneous system with the 
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Figure 2. (Colour online) Map of the local order parameter Aj/£ (upper left panel), 
the dependence of the correlation functions Ct t h (upper right panel), the STM spectra 
taken along vertical line at X—53 (lower left) and X=67 (lower right panel). The size 
of the system is 71 x 77 and there are 16% of impurities (marked by circles) which 
change the boson levels Ef locally from the homogeneous value 0.58i to Qt. The other 
parameters of this d-wave superconductor are g=0.5t, Vl mp =0, and the total number 
of carriers n=1.3. 




-0.5 



Figure 3. (Colour online) The same dependencies as in Figure [5] but LDOS is 
presented for X=17,43. We consider here the extended impurities which influence 
the positions of Ef levels of neighboring sites changing them from the homogeneous 
value 0.58£ to 0£ at the impurity site, 0.29£ at its nearest neighbor site, 0.425£ and 0.5£ 
at still further sites. 
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above parameters is superconducting. These additional impurities increase its T c . They 
also change the values of the local gaps, which are visible in the energy dependence of 
the local densities of states. 

The correlation between the position of the impurities and amplitudes of the local 
order parameters Aj and (&;) in Fig. [2] is positive and short ranged. It is important to 
note slight differences between the two correlations which are related to strictly local 
nature of the bosonic order parameter (bi) and a longer range of the fermionic one. The 
summation of Ay over the sites j neighboring i to get Aj makes the two local order 
parameters differ. The point-like character of effective interaction is responsible for the 
non monotonous dependence of these correlations on the distance. It is smooth for the 
more realistic case of longer range impurities shown in Figures |3] and HI 

The correlations between the local densities of states (at highest values) and the 
values of Aj are negative. These quantities are anti-correlated. This is true for both 
local, Figure [21 and extended Fig. [3] impurities. The correlation virtually vanishes (Fig. 
HJ) for the same kind of extended impurities introduced into non-superconducting system 
which is here modelled by the position of bosonic levels very distant E B = At from the 
Fermi energy (cf. upper left panel in figure (JTJ). 

The anti-correlation is easy to understand. Close examination of the data indicates 
that each time LDOS is measured for a site where the order parameter is large, the 
coherence peaks are relatively broad, whereas they are sharp and high in the places 
with very low or even zero value of the order parameter. For the interpretation of 
the experiments this finding means that each time the STM tip is scanning the region 
of low (or even zero) order parameter it finds sharper edges of the spectra and lower 
apparent gap (defined as the peak to peak distance in the local density of states). On 
the contrary, the spectra calculated for the patches of the sample with large values of 
the order parameters display less sharp coherence features. 

These sharp coherence-like features observed at the sites with low values of 
"effective" interactions I? resemble those found in the nonsuperconducting region 
placed inside the superconductor [63J and induced by the proximity effect. We believe 
that the same general mechanism of proximity is operating here on the local scale and 
between regions with different values of the order parameters. 

Another interesting feature of the STM spectra is their relative homogeneity at low 
energies. This qualitatively agrees with the STM data [TJ |2j [3j HJ |5] and theoretically 
a similar behaviour is also observed in the spectra obtained for the two-orbital model 
[M] with local interactions (negative Hubbard U). In this context, we have to remark 
that as mentioned earlier, calculating LDOS we replaced the Dirac delta function by 
the Lorentzian of width 77. In all calculations reported here we have taken rj = 0.04t. 
This necessarily slightly smears spectra at low energies and cuts the singularities at the 
gap edges. 

The comparison of the correlations C/^ observed in Figures [21 [3] and H] with those 
measured experimentally [TJ |2j [31 HJ |5] clearly indicates the validity of the scenario in 
which doping of the non-superconducting medium introduces bosonic impurities in a 
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Figure 4. (Colour online) Effect of extended impurities on the non-superconducting 
system with Ef=4t. The same dependencies as in Figure |3] except LDOS are shown 
for X=26,33 and the lower panels show the maps of carrier concentration: bosons 
(left), and fermions (right). The total carrier concentration is kept constant n=1.3. 
The size of the system is 71 x 77 and there are 16% of impurities (marked by circles 
on the map of Ai/t) which change the boson levels Ef from the homogeneous value 
At to Ot at the impurity site, 0.29t at its nearest neighbor site, 0.425i and 0.5i at still 
further sites. 



close vicinity of the Fermi level. 

In the above calculations the total concentration of carriers is fixed and equals 1.3. 
However, due to disorder the local values of boson n^ and fermion nf numbers vary from 
point to point. The lower panels in Figure H] show colour coded maps of nb (left panel) 
and nf (right panel). These data as well as similar maps for other cases with the fixed 
carrier concentration show that fluctuations of bosons are much stronger than those of 
fermions. In Figured] the fermion concentration varies between 0.82 to 0.88, whereas 
Ub changes from to 0.5. The relative fluctuations depend on the kind of impurities in 
the system, but as a general rule fermions fluctuate more weakly than bosons. In most 
cases studied and in agreement with the behaviour observed in Figure @] the maps of 
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Figure 5. (Colour online) Map of the local gap Aj/t (upper left panel), the dependence 
of the correlation functions Cf t h (upper right panel), the STM spectra taken along the 
vertical line at X=29 (lower left panel) and X=39 (lower right panel). The size of the 
system is 71 x 77 and there are 16% of impurities (marked by circles on the map of 
Aj/t) which change the boson levels Ef locally from the homogeneous value 0.58i to 
1.16i. The other parameters of the d-wave superconductor are g=0.5t, V* mp =0, and 
the total number of carriers n=1.3. 



concentrations show a large degree of correlations with the maps of order parameter. 

The results obtained for the second scenario according to which we introduce into 
the system bosonic impurities suppressing superconductivity are illustrated in Figure 
|5j It refers to point-like impurities, which change the value of the Ef level only 
at the impurity site. On the map of local gap one notices the diminishing of the 
superconducting gap at the impurity sites and in their close vicinity. The correlations 
between the gap value and the position of impurities take on negative values. LDOS as 
the function of u, shown in the lower panels for two cuts along X=29 (left) and X=39 
(right) show smaller fluctuations of the gap than those observed for the previous scenario 
with bosonic impurities supporting the superconductivity. The appreciable asymmetry 
of N(i, to) is due to the fact that for the parameters chosen the Fermi level of the impure 
systems has moved to the close vicinity of the Van Hove singularity. 

This scenario is not consistent with experimental findings as it leads to different 
than observed correlations between local gap and position of impurities. 



3.2. Random scatterers 

As mentioned earlier we shall assume here that impurities supporting superconductivity 
by moving the boson levels closer to the Fermi level at the same time induce potential 
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Figure 6. (Colour online) Map of the local order parameter Aj/4 (upper left panel), 
the dependence of the correlation functions Cth (upper right panel), the STM spectra 
taken along the vertical line at X=2i (lower left panel) and X=61 (lower right panel). 
The impurities are extended and influence the positions of Ef levels changing them 
from the homogeneous value 0.584 to 04 in the impurity position, 0.294 at its nearest 
neighbor site, 0.4254 and 0.54 at still further sites. In a similar way they also modify 
local energies from V i imp =0 to V? mp =— 0.29 at the impurity site, —0.1454 at its nearest 
neighbor site, —0.07754 and —0.044 at still further sites. 



scattering. Obviously this is the most realistic scenario, as in any real system with weak 
screening properties one expects charged impurities to influence not only single particle 
properties but also local interactions. 

The results are shown in Figure [6] for a system which, when undoped, is 
characterized by E B = 0.58t and low superconducting transition temperature. The 
impurities are extended and modify the bosonic levels at the central and neighboring 
sites moving them from homogeneous value 0.58t to Ot at the impurity position, 0.29t 
at its nearest neighbor site, 0.425t and 0.5t at still further sites, as in Fig. [3j We assume 
that they also modify the local "atomic" energies from y™ riv = fj to Y^ mv = —0.29 at the 
impurity, -0.145t at its nearest neighbor site, -0.0775t and -0.04t at still further sites. 
As a result LDOS show strong variations from point to point as can be seen in the lower 
panels of Figure El 

The correlation C/^ between the positions of impurities and the local values of 
the gap is positive and around 0.4 for a small distance. Again it is lower than the 
correlation of the bosonic order parameter due to the on-site nature of the latter. LDOS 
shows strong variations from site to site which are connected with the combined effect 
of changes in the local values of the effective interactions and potential scattering. It 
is the potential scattering which suppresses sharp coherence features in the LDOS. The 
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general observation that the peaks in LDOS are sharper at the sites with lower value of 
the gap remains still valid. 

The map of local values of gaps shows characteristic "two shoulder" distribution 
with the maximum of smaller gaps centered around 0.6t and for larger gaps around 
0.9t. This feature, however, is not universal and changes with the configuration and the 
extent of impurities. 

4. Summary and conclusions 

We have studied the two-component boson-fermion model in real space paying a special 
attention to the correlations between the local properties of the system as measured 
in the scanning tunnelling experiments. We have allowed the site fermion energies and 
bosonic levels to vary in a random fashion. Most attention has been paid to variations 
of the bosonic levels Ef . We propose that they are mainly affected in the process of 
doping of HTS parent compounds and this is a main source of inhomogeneity observed 
in the STM spectra. 

The presence of disorder in the bosonic subsystem provides natural explanation for 
a number of observations on HTS. In particular, the local values of order parameters are 
positively correlated with the positions of impurities (e.g. Figure [3]). Our calculations 
for short ranged bosonic impurities indicate that the sharp BCS like features observed 
in the STM spectra might not be signatures of the well developed superconducting gap, 
but rather present proximity induced structures, which appear in the regions with the 
suppressed superconducting interactions. 

In summary, our study shows the ability of the two-component boson-fermion 
model to describe the real space STM measurements (dependence of the differential 
conductance on energy, maps of the local gaps and the correlations between various 
parameters) with a high degree of accuracy. We limited ourselves here to the case of 
hard core infinitely heavy bosons. Including kinetic energy of bosons, and taking into 
account next nearest-neighbor fermion hopping will make the model more realistic and 
the results will be presented elsewhere. The model with mobile bosons has been studied 
in the homogeneous systems by e.g. [65], [66J, E7J, EH] and many other authors. This 
aspect seem to be important for the correct description of the phase diagram of HTS. 

In this paper we focused on two - dimensional lattice and mean field theory of 
superconductivity. Accordingly the questions about the dimensionality and quantum 
fluctuation around mean field solution may arise. We consider the surface layer of the 
superconductor, which is coupled to the bulk. It is easy to incorporate this coupling, 
as well as hopping to more distant sites into our approach, but this would only slightly 
change the numerical results but substantially limit the size of the clusters studied. 

The more important issue of quantum fluctuations [69] can in principle be addressed 
in a way similar to that applied recently to the inhomogeneous t-J model [70]. In the 
context of negative U - Hubbard model the inhomogeneities have been studied [71] by 
both BdG and Monte Carlo methods. The conclusion of the authors [71] is that the 
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results of the mean - field BdG calculations agree with Monte Carlo approach which takes 
into account some of the fluctuations neglected by the former method. Thus we expect 
that our results remain qualitatively correct. Influence of the phase fluctuations for the 
homogeneous version of the two-component (boson-fermion) model has been addressed 
by some of us [UJ H2] indicating the suppression of the critical temperature due to 
phase fluctuations. The energy gap (or "pseudogap") is preserved up to T*, which 
in the studied model roughly corresponds to T c obtained from the mean field BdG 
calculations. 
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